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Abstract. 

We review recent progress in analysing wave scattering in systems with both 
intrinsic chaos and/or disorder and internal losses, when the scattering matrix is 
no longer unitary. By mapping the problem onto a nonlinear supersymmetric a— 
model, we are able to derive closed form analytic expressions for the distribution of 
reflection probability in a generic disordered system. One of the most important 
properties resulting from such an analysis is statistical independence between 
the phase and the modulus of the reflection amplitude in every perfectly open 
channel. The developed theory has far-reaching consequences for many quantities 
of interest, including local Green functions and time delays. In particular, wc point 
out the role played by absorption as a sensitive indicator of mechanisms behind 
the Anderson localisation transition. We also provide a random-matrix based 
analysis of 5-matrix and impedance correlations for various symmetry classes 
as well as the distribution of transmitted power for systems with broken time- 
reversal invariance, completing previous works on the subject. The results can be 
applied, in particular, to the experimentally accessible impedance and reflection 
in a microwave or an ultrasonic cavity attached to a system of antennas. 
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1. Introduction 

Propagation of electromagnetic or ultrasonic waves in billiards 1 , compound- nucleus 
reactions [2J , scattering of light in random media and transport of electrons through 
quantum dots 0] share at least one feature in common: in all these situations one 
deals with an open wave-chaotic system studied by means of a scattering experiment, 
see figure 1 for an illustration. Here, we have a typical transport problem where 
the fundamental object of interest is the scattering matrix S, which relates linearly 
the amplitudes of incoming and outgoing fluxes. However, under real laboratory 
conditions there is a number of different sources which cause that a part of the 
flux gets irreversibly lost or dissolved in the environment. As a result, we encounter 
absorption and have to handle the iS-matrix, which is no longer unitary. Statistics 
of different scattering observables in the presence of absorption are nowadays under 
intensive experimental and theoretical investigations, starting from early experiments 
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Figure 1. A sketch of a typical experimental setup with microwave billiards. A 
flat chaotic cavity is feeded with microwaves through an attached coaxial cable 
(i.e. a scattering channel). On average, 1 — T part of the incoming flux, where 
T < 1 is the so-called transmission coefficient, is reflected back directly from the 
cable-cavity interface (port) without exciting long-lived resonances in the cavity. 
If the cavity is thin enough then only a transverse electric wave can propagate 
inside. The electric field has only a vertical component, which is uniform in 
vertical direction and distributed nontrivially in the plane. Therefore, there is 
a voltage between plates as well as a current due to the in-plane magnetic field. 
The impedance is a quantity which relates linearly the port voltage to the port 
current. Fluctuations of eigenmodes and eigenfrequencies result in fluctuations of 
the impedance or S-matrix, as the driving frequency or port position is changed. 



on reflection and energy correlations of the S- matrix [5] |S] . More recently, total cross- 
sections [7], distributions of reflection [5] and transmission |5] coefficients as well as 
that of the complete S matrix |1(J | in microwave cavities, properties of resonance 
widths 1 1 1 1 in such systems at room temperatures, dissipation of ultrasonic energy 
in elastodynamic billiards and fluctuations in microwave networks became 
experimentally available. Theoretically, statistics of reflection, delay times and related 
quantities were considered first in the strong J3| an d then weak ^3] absorption limits 
at perfect coupling and very recently at arbitrary absorption and coupling for several 
symmetry classes dU H3 HHJ dl E01 EU E2 E3 EU 

Another insight to the same problem comes from looking at it not from the 
"outside" , but rather from the "inside" . Then the prime object of interest turns 
out to be the impedance Z relating linearly voltages to currents at the system 
input , see figure 1 . By proper taking into account the wave nature of the 

current fZ7\ |2H| the cavity impedance can be seen as an electromagnetic analogue of 
Wigner's reaction matrix of the scattering theory. This can be easily understood 
qualitatively through the well-known equivalence of the two-dimensional Maxwell 
equations to the Schrodinger equation, the role of the wave function being played 
by the electromagnetic field (the voltage in our case). Then the definition of the 
impedance becomes formally similar to the definition of the reaction matrix (which 
relates linearly the scattering wave function to its normal derivative on the boundary). 
The impedance is, therefore, related to the local Green's function of the closed cavity 
and fluctuates strongly due to chaotic internal dynamics. 

The imaginary part of the local Green's function (which is proportional to the real 
part of Z) is well known as the local density of states (LDoS) and has a long story of 
studies in disordered electronic systems, see for a recent review. Actually, a closely 
related quantity emerges in the context of spectra of complex atoms and molecules 
where it has the meaning of the total cross-section of indirect photoabsorption [23 (see 
also [SUES)- It also appears in studies on spontaneous light emission by atoms placed 
in chaotic cavities [SB.. As to the real part of the Green's function, it seems to have 
no direct physical meaning in mesoscopics while it has the meaning of reactance in 
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electromagnetics, where both real and imaginary parts are experimentally accessible. 

In this review we discuss an approach developed recently by us in short 
communications [211 1221 &\ which treats both inside and outside aspects of the 
problem on equal footing. In this capacity it provides a uniform and deeper 
understanding of various results on absorptive scattering obtained earlier in ^3 El 
I18| . In particular, the method allows one to study very efficiently the distribution of 
the local Green's function (complex impedance) at arbitrary absorption and to relate 
it to that of reflection, thus linking somewhat complementary experiments and 
|2"5) together. Although calculations are most explicit and simple for fully chaotic 
systems, when one can rely upon the random matrix theory (RMT), our method 
actually has relevance in a much broader context beyond RMT that involves many 
interesting aspects of disordered mesoscopic systems with absorption, including effects 
of the Anderson localisation. From that point of view, the method opens an attractive 
possibility to look at some long-standing problems (e.g. statistics of time delays) from 
a different perspective, see [35] as well as subsections 4.1.6 and 4.1.7 below. 

2. Reflection, time delays and resonance spectrum 

In this section, we provide a short description of the scattering approach to the 
problem. The resonance energy dependence of observables becomes explicit in the 
well-known Hamiltonian approach to quantum scattering, which was developed first 
in the context of nuclear physics [21 EH1 SHI an d can be easily adopted for models 
emerging in quantum chaotic scattering and mesoscopic physics, see e.g. [H I3U1 HU1 14T] 
for reviews. This framework turns out to be also most suited to take a finite absorption 
into account. The starting point is the following fundamental relation between the 
resonance part of the scattering matrix and the Wigner's reaction matrix K: 

s{e) = \ + Tk{e) ' k{e) = \v\e-h)- 1 v, (1) 

The Hermitian Hamiltonian H of the closed system gives rise to N real energy levels 
(eigenfrequencies). Those are coupled to M continuum channels via the NxM matrix 
V of coupling amplitudes V£ (n = 1 . . .N, c = 1 . . . M), and as a result are converted 
to N complex resonances. To see this, we expand Q in a Taylor series in K and, 
after regrouping the terms, bring the resulting expression to another well known form 

S{E) = l-W \ 1 V, W eff = H - IVV^ (2) 

for the S-matrix. The effective Hamiltonian 7Y c ff emerging here characterises the 
open system and is the non-Hermitian counterpart of H . The factorized structure 
of the anti-Hermitian part is necessary to ensure the unitarity of S. The coupling 
amplitudes V change very slowly with the energy (far from the channels thresholds) 
and one can safely consider them to be energy-independent. In such a resonance 
approximation the complex eigenvalues E n = E n — ^T n of 7i c ff, with the energies E n 
and escape widths T„ > 0, are the only singularities of the S matrix in the complex 
energy plane. As required by causality 02 > they are located in the lower half plane 
and correspond to the long-lived resonance states formed on the intermediate stage 
of a scattering process. The corresponding (left and right) eigenvectors form the so- 
called bi-orthogonal system. Recent discussion concerning applicability of the effective 
Hamiltonian approach to potential scattering problems (like those with cavities) can 
be found in g3 |H gSJ HH . 
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Flux conservation requires S(E) to be unitary at the real values of E. It is useful 
to define at real f2 the following matrix 

Rn(E) =Sfi(E- $Sl)S(E + = 1 + iffl^ 1 1 V . (3) 

E — 5S2 — W^j £ + 2" _ "off 

The second equality here results from the substitution of @ and then making use of 
the identity iVV^ + tt = (E+ ± Q - W c fr ) - - §0 - ) which leads to cancellations 
of the cross-terms. Expression (01 tends to unity as f2 — ► and it is the Wigner-Smith 
time-delay matrix gaEHj Q{E) = -\S\E)^ E S(E) = -i^R n (E)\ n=0 (we put H = 1) 
which determines the unitarity deficit of Rq to the linear order |49| : 

R n (E)^l + iQQ(E)+0(n 2 ), Q(E) = V \ 1 \ V. (4) 

— rtcffj 1 & — ties 

Such a factorized representation |5()| for the time-delay matrix (which contains 
no longer the energy derivative) is a consequence of the resonance approximation 
considered. It serves to make a connection of time delays to the resonance spectrum 
most explicit. The matrix element Q cc — (b^b) cc may be physically interpreted as 
the scalar product (or "overlap") between the internal parts b c — (E — H. c tf)~ 1 V c of 
the scattering wave functions for waves incident in the channels c and c', respectively 
|5U| . In particular, the mean time delay in the channel c given by the diagonal element 
Q cc coincides in such an approximation with the dwell time given by the norm of b c . 
One should distinguish generally between Q cc and the so-called proper time-delays 
q c (eigenvalues of Q). Taking their sum, one comes to a weighted mean time delay 
characteristic Q w = jftrQ — jjJ2 c Q cc = i?Sc9c, the so-called Wigner's time 
delay which is known, see e.g. [220^10^1 to be determined by the energy derivative 
of the total scattering phaseshift: Q w = — jj -Bs log det S. Diverse aspects of delay 
times in quantum chaotic scattering |40| as well as in a general quantum mechanical 
context [521153] can be found in the cited literature and references therein. 

As is well known, statistics of spectra of closed quantum systems with chaotic 
classical counterparts are to a large extent universal and independent of their 
microscopic nature. This remarkable universality provides one with a possibility to 
use the RMT ,54 for an adequate description of many physical properties of such 
systems |55j . According to the general paradigm we replace the actual Hermitian 
part H of the effective Hamiltonian J5J) with a random Hermitian matrix H taken 
from one of the three canonical Wigner-Dyson's ensembles labelled by the symmetry 
index [3 according to the symmetry of the system under consideration. The Gaussian 
orthogonal (GOE, (3—1 and H real symmetric) and unitary (GUE, = 2 and H 
Hermitian) ensembles stand for systems with preserved or fully broken time-reversal 
symmetry (TRS), respectively. The remaining Gaussian symplectic ensemble (GSE, 
(3 = 4 and H self-dual quaternion) is relevant for description of time-reversal systems 
with strong spin-orbit scattering. The limit of large N — > 00 is supposed to be finally 
taken. Then eigenvalues are distributed on the finite interval according to Wigner's 
semicircle law, which determines locally the mean level spacing A. The most appealing 
feature of the RMT approach is that quantities related to spectral fluctuations when 
expressed in units of A ("unfolding") do not depend on microscopic details (i.e. the 
particular form of the distribution of H or the profile of A) and become uniform 
throughout the whole spectrum |55| . For practical reasons we thus restrict ourselves 
to considering fluctuations at the center of the spectrum (E — 0) only. Similarly, 
the results turn out to be also independent of particular statistical assumptions on 
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coupling amplitudes V£ as long as M C JV |561 H§| . The amplitudes may be chosen 
as fixed [2] or random |38) variables and enter final expressions only in combinations 
known as transmission coefficients (also sometimes called sticking probabilities) 

T c =i-m> = , Kc = ^p ( 5 ) 

1 1 (l + Kc) 2 2NA v ; 

where S cc stands for the average (optical) S matrix. The transmission coefficients 
are assumed to be input parameters of the theory, the cases T c <§; 1 and T c = 1 
corresponding to an almost closed or perfectly open channel "c" , respectively. 

Absorption is usually seen as a dissipation process, which evolves exponentially 
in time. Strictly speaking, different spectral components of the field may have 
different dissipation rates. However, frequently this rather weak energy dependence 
can easily be neglected as long as local fluctuations on much finer energy scale ~ A are 
considered. As a result, all the resonances acquire one and the same absorption width 
r > additionally to their escape widths r„. The dimcnsionless phenomenological 
parameter 7 = 27iT/A characterizes then the absorption strength, with 7 < 1 or 7 > 
1 corresponding to the weak or strong absorption limit, respectively. Microscopically, 
it can be modelled by means of a huge number of weakly open parasitic channels |rjll57| 
or by coupling to a very complicated background with almost continuous spectrum 
[lfij . see also |58| . In microwave billiards such an approximation is frequently very 
good to account for uniform Ohmic losses which happen everywhere in non-perfectly 
conducting walls. However, in some experimentally relevant situations as, e.g., 
complex reverberant structures |59| or even microwave cavities at room temperature 
|11II46| an approximation of uniform absorption may break down, and one should take 
into account instead localized-in-space losses which will result in different broadenings 
of different modes. The latter are easily incorporated in the model by treating them 
as if induced by additional scattering channels, see e.g. [T2|. An alternative scheme 
of treating localised-in-space surface absorption is discussed in [221 • (Discussion of a 
formal theory of scattering for complex absorbing potentials can be found in |f>()|.) 

Operationally, the uniform absorption can equivalently be taken into account by 
a purely imaginary shift of the scattering energy E — > E + j^Y = E 11 so that the 
5-matrix S 7 (E) = S(E 1 ) becomes subunitary. The reflection matrix R 1 = S^S^ 
provides then a natural measure of the mismatch between incoming and outgoing 
fluxes. It can be obtained from i?n, by analytic continuation in O from a real to 
the purely imaginary value fl = iY, yielding [16| : 

Ry(E) = Rn=ir(E) = 1 - YQ 7 (E) , Q (g) = yt 1 1 y. ( 6 ) 

(ii 7 — rtcff) 1 Jii-y — rtcs 

This representation is valid at arbitrary value of Y. In the limit of small Y one can 
neglect the difference between Q 7 and Q, resulting in the approximate expression 
j^E] i? 7 — 1 — TQ following from 10} . It is therefore tempting to keep for Q 1 the 
meaning of the time-delay matrix at finite absorption as well (see, however, discussion 
in By construction, i? 7 is a Hermitian matrix, its positive reflection eigenvalues 

r c = 1 — Yq c < 1 are related to proper time-delays q c (at finite absorption) . 

Considering quantum scattering with no internal dissipation, Y = 0, the ^-matrix 
unitarity ensures that the reflection coefficient R c — \S CC \ 2 in any given channel is 
simply related to the quantum mechanical probability to exit via any of the remaining 
channels, known as the transmission probability: R c = 1 — J2b^c \S bc \ 2 - For an 
absorptive system the last equality is violated, but still the quantity r c = 1 — R c 
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can be interpreted as the quantum mechanical probability that a particle entering 
via a given channel never exits through the same channel. Hence it was suggested 
to call r c the probability of no return (PNR) |17j . In the particular case of a single 
open channel fluctuations of the PNR r arise solely due to absorption (neglecting 
dissipation trivially results in r = 0). At weak absorption r w Tq , so that the PNR 
is just simply proportional to the time-delay. 

At last but not least, the matrix Z = iK(Ej) has the meaning of the normalized 
cavity impedance in such a setting, see [23 EHj for further details. 

3. Correlation functions 

Any observable in chaotic resonance scattering exhibits strong fluctuations over a 
smooth background as the scattering energy or other external parameters are varied. 
Usually these variations occur on two essentially different energy scales. This fact is 
conventionally taken into account by decomposing fluctuating quantities into a mean 
part and a fluctuating part, the former being understood as the result of spectral 
or (assumed to be equivalent) ensemble average (•••). In this section we consider 
statistics as determined by a two-point correlation function of the fluctuating parts 
(frequently called a "connected" correlation function) : ( AB) colm = (AB) — (A) (B) . 
We restrict ourselves below to the cases of preserved (/3 = 1) and broken ((3 = 2) TRS 
(the symplectic case (3 = 4 proceeds along the same lines). 

3.1. Impedance 

Let us start with considering the simplest case of the impedance correlations. The 
problem can be fully reduced to that of spectral correlations determined by the two- 
point cluster function Y 2 ^(uj) = S(uj) — A 2 {p{E{) p(E 2 )) conn , where to — (E 2 — Ei)/A 
and p(E) being the level density. It is easy to satisfy oneself that for the mean 
impedance at E = holds (Z ab ) = n a 5 ab . To calculate the energy correlation 
function Cf cd (cu) = (Z ab * (E 1 )Z cd (E 2 )) conn , it is instructive to write Z ab (E) in the 
eigenbasis of the closed system: Z ab (E) = ^J2n v n* v n/( E - E n + %T). A rotation 
that diagonalizes the (random) Hamiltonian matrix H transforms the (fixed) coupling 
amplitudes V£ to Gaussian distributed random coupling amplitudes nJJ with zero mean 
and covariances (wjj'wj,) = (2K a A/ir)S ab 5 nm . In such a representation the energy 
correlation function acquires the following form: 

c abcd {uj) = y\ ( — - — , - — — ) (7) 

and the averaging over coupling amplitudes and that over the spectrum can be done 
independently. The Gaussian statistics of v results in 

1 / 7r \2 / a b* c* d \ cab red , / cac cbd , c cad cbc\ r /o\ 

4(3-) {v n v n v m v m ) = K a n c d d + K a K b (d d + difjd d )d nm (8) 

where 5\p term accounts for the presence of TRS, when all are real and Z is symmet- 
ric. It is useful to represent the spectral correlation function in a form of the Fourier 
integral /" dt, / °° dt 2 e- r ^ +t ^/ 2 e iE(t ^ e^-E^+t,)^ ^i(E„i 1 -S m t 2 )^^_ Due 
to the uniformity of local fluctuations in the bulk of the spectrum, one can integrate 
additionally over the position E of the mean energy: J -^-e 1 ^* 2- * 1 ) = -^■<5( * 2 t ~ tl ), 
where tu = 2n/A is the Heisenberg time. It is natural therefore to measure the 
time in units of tn- From the known RMT spectral fluctuations one also has for 
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n^m (1 - N) (e 27ri(£ "~ Bm)t/A ) conn = b 2 ^(t), where b 2 ,p(t) is the spectral form factor 
defined through the Fourier transform of Y 2 ^(u>) — dte 2muJt b2,B(t) 

&2,/3=i (t) = [1 - It + 1 log(l + 2i)]0(l - *) + [t log(i±i) - l]0(t - 1) (9a) 

6a,/3=2(t) = (l-t)Q(l-*) ( 96 ) 
at t > and b 2 ,g(—t) = b 2 ,g(t). Combining all these together, we arrive finally at 

poo 

Cf cd (uj) = / dte 2niu)t Cf cd (t) (10a) 
Jo 

Cf cd {t) = 4 e-i" [« « c [l - 6 2)/3 (t)]5 ob <r d + K a /^ oc <5 bd + 5 1/3( 5 ad ,5 bc )] . (106) 

Similar in spirit calculations were done in a context of reverberation in complex 
structures in jB2 and in a context of chaotic photodissociation in 63, (pj. 

The form factor (|10foJ) is simply related to that of K matrix elements at zero 
absorption as C^ bcd (t) = e~ l1 C^ cd {t) . Such a relationship between the corresponding 
form factors with and without absorption is generally valid for any correlation function 
which may be reduced to the two-point correlation function of resolvents (see \7\ and 
the discussion below, e.g., for the case of the S matrix). This can be easily understood 
as a result of the analytic continuation 2itlo — > 2ttuj + ij of the energy difference u 
reflecting switching on the absorption (see the previous section). 

The obtained expressions describe a gradual loss of correlations in Z matrix 
elements as the energy difference grows; generally, Cz-(oj— >oo) — > 0. At uj = 0, ((lOajl 
provides us with impedance variances Cf ab (0) = var(Z ab ) = (\Z ab \ 2 ) - \(Z ab ) | 2 , which 
were recently studied experimentally in [HS]. In analogy with the so-called elastic 
enhancement factor considered frequently in nuclear physics |66) . one can define the 
following ratio of variances in reflection (a=b) to that in transmission (a^b): 



Wz,B = ^ ' V >- = 2 + 8 w - / dse- s b 2 ^) (11) 



o 



v /var(Z aa )var(Z bb ) 
var(Z afc ) 

where the second equality follows easily from l|106t (note that the coupling constants 
n a .b are mutually cancelled here). Making use of &2„s(oo)=0 and b 2 ^(0)=l, one can 
readily find Wz,b in the limiting cases of weak or strong absorption as: 

f 2 + 5x8 at 7 < 1 

Wza = < ■ 12 

\l + S 10 at 7 >l V ' 

Wz,b decays monotonically as absorption grows. In the case of unitary symmetry, 
<|96(1 and l|ll|) yield explicitly Wz,2 = 1 + -(1 — e r ) in agreement with [SS]. It is 
hardly possible to get a simple explicit expression at finite 7 in the case of orthogonal 
symmetry. However, a reasonable approximation can be found if one notices that 
the integration in 1)11(1 is determined mainly by the region s < 1, so that one can 
approximate &2,i(s) ~ (1 — 2s + 2s 2 )0(l — s) through its Taylor expansion. Performing 
the integration, one arrives at Wz,i ~ 3 — 7~ 2 [(4 + 7 2 )(1 — e~ 7 ) — 27(1 + e -7 )], 
which turns out to be a good approximation to the exact answer at moderately strong 
absorption (deviations are seen numerically only at 7 ~ 1). 

3.2. Scattering matrix 

The energy correlation function of the scattering matrix elements 

/>oo 

Cf cd (uj) = (Sf*(E 1 )S c 7 d (E 2 )) conn = / dte^Cf^it) (13) 

Jo 
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is a much more complicated object for an analytical consideration as compared to 
(7J . The reason becomes clear if one considers again the pole representation of the S 
matrix which follows from ©: S ab (E) = S ab - iJ2 w^w^/ (E - £ n ). Due to unitarity 
constraints imposed on S at real E the residues and complex resonance energies 
develop nontrivial correlations |2Hj ■ Although results on statistics of resonances in the 
complex plane @H1 EH EH EH] as well as the corresponding residues EHl EH became 
recently available in some particular cases, this knowledge is insufficient as yet for 
calculating S'-matrix correlations in general. The separation like (J7J into a "coupling" 
and "spectral" average is no longer possible and can be done only by involving some 
approximations |71|. The powerful supersymmetry method |721 1^] turns out to be 
an appropriate technique to perform the statistical average in this case. In their 
seminal paper [2], Verbaarschot, Weidenmuller and Zirnbauer performed the exact 
calculation of l|13|) at arbitrary transmission coefficients (and zero absorption) in the 
case of orthogonal symmetry. This finding was later adopted [7] to include absorption. 
The corresponding exact result for unitary symmetry has been recently presented by 
us in [32] (see also [73] concerning the S'-matrix variance in the GOE-GUE crossover 
at perfect coupling) and is discussed below. 

The calculation proceeds along the same lines as in [2]. The final expression 
for both the connected correlation function and its form factor can be equally 
represented as follows: 



Cf cd = 6 ab 6 cd T a T c ^/(l - T a )(l - T c )J ac + (S ac S bd + S w 6 ad S bc )T a T b P ab . (14) 

Here, the 8\p term accounts trivially for the symmetry property S ab — S ba in 
the presence of TRS. J ac and P ab defined below are some functions (of the energy 
difference lo or the time i), which depend also on TRS, coupling and absorption 
but already in a nontrivial way. As a result, the elastic enhancement factor 
Ws t /3 = \J var(S aa )var(S bf, )/var(S ab ) is generally a complicated function of all these 
parameters, in contrast to In the particular case of perfect coupling, all T c = 1, 

one has obviously from l|14|) that Ws,p = 1 + <5i/3 at any absorption strength. 

We consider first expression (|14fl in the energy domain at real uo (no absorption) . 
The functions J ac (uj) and P ab (w) can be generally written as certain expectation 
values in the field theory (nonlinear zero-dimensional supersymmetric cr-model) whose 
explicit representations depend on the symmetry case considered, we refer the reader 
to [7J] for general discussion. In the (3 — 1 case of orthogonal symmetry the well- 
known result of [2] reads as follows: 

p ( \ _ 1 1 MH-Mi) , M2(l+A*2) i Mo(l-Mo) \ -r \ r h \ 

r ab \u) - ^ (1+TaAIi)(1+T(;Mi) -t- (i+T aAI2 )(i+T bM2 ) + (i-T a/ i )(i-r !)A (o)K M / A( y 100 ) 

with T M = n e [ (i+T 1 /jO(i+T being tne so-called "channel factor", which 

accounts for system openness, and ((• • is to be understood explicitly as 

a °W d/x 2 / d/ioTT— - — — 1/2 — - — rj- — - — rj (...)• (16 

8 Jo Jo Jo [(1 + M1JM 1 + V2)M / [Vo +Mi) (Mo + 

In the [3 = 2 case of unitary symmetry we have found |34j that 

■M") = ((rfc + rfc) (tt%t + rfc)^ (17a) 
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1 



d Ml / d/i (/ii + W>r 2 e [2 ™^+^ (...)• (18) 



o 

In the important particular case of the single open channel (elastic scattering), the 
general expression for the (3 = 2 case simplifies further to 

(rwsm)^ - ^ /°V / V! ,, \t+t,,T ^ { ^ o] • (i9) 

Jo Jo Mi + Mo U + J Mir 
Finally, putting above lu — ► lu + i"f /2ir accounts for the finite absorption strength 7. 

To consider <|14|) in the time domain, i.e. the form factor Cg bcd (t), we notice 
that the variable t = 5(^1 + M2 + 2/io) f° r /? = 1 or t = /ii + /io for /3 = 2 plays 
the role of the dimensionless time (in units of tn). The corresponding expressions for 
P a b(t) and J ac (t) can be investigated using the methods developed in [Sj3 EH EI] For 
orthogonal symmetry it was done in [7J, where the overall decaying factor e~ 7t due to 
absorption was also confirmed by comparison to the experimental result for the form 
factor measured in microwave cavities. It is useful for the qualitative description to 
note that P a b(t) and 2J ac (t) are quite similar to the "norm leakage" decay function [75] 
and the form factor of the Wigner's time delays 02] > respectively (they would coincide 
exactly at 7 = 0, if we put T a j, c = appearing explicitly in denominators above). 
Then one can follow analysis performed there, see also to find P a b(t) « e -7 ' and 
Jac(t) ~ {2t/ (3)e~ lt as exact asymptotic at small times [7T], while they both become 
proportional to e~ J H~ M ^/ 2 ~ 2 at large times. 

Such a power law is characteristic for open systems [35| E3 ■ Physically, it 
results from width fluctuations, which diminish as the number M of open channels 
grows In the limiting case M — > 00 and T c — ► 0, all the resonances acquire just 

the same escape width ^2 C T C (in units of tjj 1 ), which is often called the Weisskopf's 
width [75], so that the total width is 7t = Tl c Tc + 7- Then there occur further 
simplifications: P a b(t) = e~ 7T * and J a b{t) = [1 — &2. / g(^)]e _7T *, that results finally in 

Cf cd {u) = [ ° ° + ° lfi ° ' 11 b + 6 ab S cd T a T c / dt[l - b 2 M^e-^-^K^O) 
7T - 2tt\u J 

For the case of (3 = 1 this result (at zero absorption) was obtained earlier by 
Verbaarschot [HHj. In the limit considered, expression (|20|l is very similar to l|10o|) . 
IjlOfcll , so that the enhancement factor Ws.p is given by the same (|11H where 7 is to be 
substituted with jt- At 7t 2> 1 (large resonance overlapping or strong absorption, 
or both) the dominating term in H2U|) is the first one, which is known as the Hauser- 
Feshbash relation J77J, see [751 r751 f5D| for discussion. Then W s ,p = 2/f3 = W z ,p that 
can be also understood as the consequence of the gaussian statistics of S (as well as 
of Z) in the limit of strong absorption |8(J|. 

3.3. Reflection 

In this subsection we consider fluctuations in the weighted-mean reflection coefficient 
defined as r = i tri?. Its average value (r) was recently calculated in ^SJ, where the 
following exact result was found to be 



(r) = 1 - — 

V 1 M 



poo 

1-7 / dte'^Pit) 
Jo 



(21) 
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relating (r) to the so-called "norm- leakage" function P(t) introduced in [75]. The 
actual value of (r) changes with growing absorption between 1 — at 7 <C 1 and 
1 — 17 I] c T c at 7 > M. For the sake of simplicity we consider the case of M equivalent 
channels below, all T c = T. For large M ^> 1 the widths cease to fluctuate, yielding 
P{t) = e~ MTt , so that JUJ results in (r) = (MT + 7(1 - T))/(MT + 7) [HI. 

The correlation function C«(u;) = {r{Ei)r(Eh) = T 2 (Q w (E 1 )Q w (E 2 } coun is 
simply related by virtue of © to that of Wigner's time delay Q w = jj tr Q 1 (at finite 
r). Amounting to a four-point correlation function of resolvents, its evaluation is 
generally beyond the present day state of art in the supersymmetry method. However, 
by approximating Q w sa -j J2 n ^n/[(E- E n ) 2 + ^(T n +T) 2 ], cf. JBJ, and making use of 
the rescaled Breit-Wigner approximation (RBWA) of Gorin and Seligman j7J|, one can 
find an expression which is supposed [7J to work reasonably well at finite absorption. 
The final result can be naturally represented in the form of the Fourier-integral 

2t 2 f°° 

C R (u) = ^J dtcos(27Tu;t)e-^F R (t) (22) 

where the form- factor (we have set apart in (1221) the trivial absorbing factor e _7t ) is 

F (t\ = ° 2 ^ - hu3(t) Cl ^ )2 

R[ ' (1 + 2Tt//3) A W2 (1 + Tt/(3) M P 1 ' 

with the function c n ^(t), n — 1,2, defined as follows: 

/-oo n-1 ~{l3 1 /2T){l+2Tt/p)x 

CnAt) = J! + I) I d X (1 - x)M0/2+n . (24) 

It is instructive to consider the limiting cases of weak and strong absorption in 
more detail. At 7 <C T < 1, the 7-dependence and thus ^-dependence of c n) p(t) is very 
weak, so that c n ^{t) w 1 for t -C I/7. As a result, the form-factor reduces to 

F R (t) = (1 + 2Tt//3)" M ' 3 / 2 - b 2 ^(t)(l + Tt/p)- M ? (25) 

which is, as expected, just the form-factor of the Wigner's time-delay correlation 
function at zero absorption calculated in the same scheme RBWA. In the opposite 
limiting case 7>T, the x-integration in c n ^{t) is determined mainly by the region 
of small x < T/7 1, that yields readily 



f MT 
V 7 



2 r 



1 + 2/M0 b 2 , p {t) 



(26) 



(1 + 2Tt/(3) M f } ' 2 + 2 (1 + Tt/f3) M f 3 + 2 _ 

Comparing the large time asymptotic of F R (t), one sees that a crossover from £~ A W 2 
to £- M /V2-2 behavior happens at 7 ~ T, as we go from weak to strong absorption 
regime. This behaviour is indicative of an additional decay of correlations induced by 
strong absorption on top of the pure exponential one. It is most pronounced for a 
single-channel cavity with TRS, corresponding to changing from t -1 / 2 to t~ 5 / 2 . 



4. Distribution functions 

4-1- Single-channel scattering 

If.. 1.1. Reflection coefficient and the local Green function. It is quite clear, that the 
joint distribution of the real and imaginary parts of K 

V{u,v) = (5{u-ReK)5(v + lmK)) (27) 
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determines fully the statistics of the impedance and the S'-matrix. By definition 
K = k(NA/tt)Gu is related to the diagonal element Gu = (i|(-E 7 — of 
the Green function of the closed system in the position representation, taken at the 
position i of channel attachment. The RMT assumption for H implies however that 
the basis of representation may be chosen arbitrary. One studies, therefore, statistics 
of the dimensionless local Green function K = u—iv, with v > being the local density 
of states (LDoS). We consider first the case of perfect coupling n — 1 (T = 1). In 
particular, this choice ensures the normalisation of the mean LDoS such that (v) = 1. 
In this case, the following general form of the distribution function 

V(u,v) = ^P (x), x= u2+ £ + 1 >l. (28) 

may be easily established [2TJ . It initially emerged in |81l l%2] in the course of tedious 
calculations for the (3=2 symmetry class, but neither origin nor generality of such 
a form were appreciated. Actually, the validity of the form i|28[) follows directly 
from exploiting the two following fundamental properties of the S'-matrix statistics 
at perfect coupling: (i) the uniform distribution of the scattering phase 9 £ (0, 2tt); 
and (ii) the statistical independence of 9 and the S matrix modulus. In the cases of 
chaotic systems with pure symmetries, both these properties can be verified making 
use of methods of |57]. In this paper we are proving directly by an alternative, 
more powerful method [22j . It holds much beyond the universal RMT regime for a very 
broad class of disordered Hamiltonians, in particular for the case when localisation 
effects already play an important role. Moreover, this method will help us to verify 
that i|28J) holds in the crossover regime between the pure ensembles. For the time 
being the uniformity of the scattering phase distribution will be taken for granted. 

Substituting K = u — iv into Q , we can immediately infer that the variable x is 
directly related to the reflection coefficient, parameterizing it in the following way: 

x — 1 u 2 + v 2 — 1 

- , cot 9 = — 

1 2u 

The joint distribution V(u,v) of u and v and that P(x,9) of x and 9 are related 
by means of the Jacobian | gfen | = v~ 2 as can be verified by a straightforward 
calculation. Finally, P(x,9) = j^Pq(x) results from the uniform distribution of 
the phase and in this way Po(x) acquires the physical meaning of the normalized 
distribution of the reflection parameter x. 

Relationship i|28[l is one of our central results, which despite its apparent simplicity 
will have many far-reaching consequences which are not easy to guess otherwise. Let 
us note the invariance property of this distribution under the change iK — > 1/iK, 
meaning that both the impedance and its inverse (i.e. the admittance) must have one 
and the same probability distribution. Another remarkable feature of H28J) is, that it 
relates the distribution of the local Green function in the completely closed system to 
the distribution of the reflection coefficient in the perfectly open one. This fact will 
be fully utilised later on for extracting statistical properties of the Wigner time delay. 

It is instructive to discuss first the exact limiting statistics 

, (x+l)^ 2 e -«(s+l)/2 ( a< ^l) 

, , i 2r(/3/2+i) \ i ) \a<^ i) 

P {x) = <( a = |/?7, (30) 

a —a(x—l)/2 



S 7 = v^e ie , r = ^-^, cot#=^£ -. (29) 



2 



e'^- 1 " 2 (a»l) 



which can be found |21) in the weak or strong absorption limit. In the former case, one 
can follow |611 ITS] in using a perturbation theory to relate reflection to the time delay 
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in the ideal cavity without absorption (see © thereafter). The latter quantity has a 
known distribution 001113 V t (t = = [(|)/V2/r(|)] T -/3/2-2 e -/V2r , which yields 
the first line above. In the opposite case of strong absorption, the real and imaginary 
parts of S'-matrix become statistically independent Gaussian distributed variables 
|8DI IT4] . resulting in the so-called Rayleigh distribution P(r) ~ (7/3/2)e~ I " y ' 3 / 2 of 
reflection valid at r « i(ar — 1) -C 1, and the second line in l[3U)l follows. 

As to arbitrary absorption, an exact result was obtained first for unitary symmetry 
((3 = 2) by Beenakker and Brouwer [TJI- It is convenient to keep the scaled absorption 
parameter a, representing their result in the following form 

P (x) = ^7f[ A 00 + l)/2)^ /2 + s] e- a ^ +1 )/ 2 (31) 

with a-dependent constants A = e Q — 1 and _B = l + a — e a , and with A/2 = 1 standing 
for the normalization. For the case of symplectic symmetry (/3 = 4), the exact form 
was reported very recently |21| . The explicit derivation is given in Appendix A, the 
final result being: 

PT(x) = P guc (x) + [y + I) 2 - 7(7 + l)0r + 1) + 7] e"^ jfdi^p (32) 

where P,f ue (a;) i s the distribution l|31|) for /3 = 2 taken, however, at a = 27. The 
case of orthogonal symmetry, the most welcomed experimentally, turns out to be 
extremely tricky. It was suggested in [21], see also ]W\, that expression l|31[) (with 
N fj = a/ (AT (P/ 2+1, a) + Be~ a ) and T(v,a) = J^dpp^e-P) is an appropriate 
interpolation formula at (3 — 1 . It incorporates correctly both known limiting cases of 
weak or strong absorption and a reasonable agreement with available numerical and 
experimental data was found there in a broad range of the absorption strength. Below 
we provide an exact analytical treatment of this case by the method suggested in [221 • 

4-1-2. V(u,v) and the spectral correlation function of Green's function resolvents. 
We establish now the general relation [221 between the joint distribution function (|28|l 
at arbitrary finite absorption and the energy autocorrelation function 

° n ' Z+) = ( z_ - iO - K (E - p - iO) z + + iO - K (E + p + iO) ) (33) 
of the resolvents of the local Green function Kg at zero absorption (T = 0, indicated 
explicitly in this subsection with the subscript "0" ) . Distribution (|28|l can be obtained 
from l|33ll by analytic continuation in O from a real to purely imaginary value = IT 
as follows. Kq(E) is an analytic function of the energy in the upper or lower half- plane 
and, thus, can be analytically continued to the complex values: Kq(E ± ^T) = u + iv, 
v > 0. This allows us to continue then analytically the correlation function l)33|) from 
a pair of its real arguments to the complex conjugate one: z+ = (z_)* = z' + iz" , 
z" > 0. As a result, function (|33|l acquires at = ir the following form: 

C(z',z") = Cn=ir(z-,z+) = (-r-, \2~Ti « 

\ (z — u) A + (z" + vy 

V(u,v) 
(z< - uf + (z" + vf 

The second line here is due to the definition (|28|l . To solve this equation for V(u,v), 
we perform first the Fourier transform (FT) with respect to z' that leads to 

C(k, z") = J dz'e lkz C(z', z") = J dvV(k, v) z „ + v (35) 



du dv— r^TT, ro- (34) 
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where V(k,v) is the corresponding FT of V(u,v), Being derived at z" > 0, Eq. 
can be analytically continued to the whole complex z" plane with a cut along negative 
Rez". Calculating then the jump of C(k,z") on the discontinuity line z" = —v 
(v > 0), we finally get the following expression 

V(k, v) = ^$ \d(k, -v - iO) - C(k, -v + i0)l (36) 
2tt z \ L J 

with the Heaviside step function 0(i>). The inverse FT of l|36(l yields V(u, v). 

This completely general relationship is another our central result. It resembles 
(and reduces to) the well-known relation between the spectral density of states and 
the imaginary part of the corresponding resolvent operator, when the case of one real 
variable is considered. In contrast, the case of the distribution of two real variables 
requires to deal with the two-point correlation function. Physically, the latter is a 
generalized susceptibility describing a response of the system under consideration. 
This fact suggests to consider our formula in a sense as a "fluctuation dissipation" 
relation: The l.h.s. of l|36|l stands for the distribution (of K) in the presence of 
dissipation / absorption whereas the correlation function (of resolvents of K) in the 
r.h.s. accounts for fluctuations in the system, i.e. for arbitrary order correlations in 
the absence of absorption. 

The main advantage of the derived relation is that the correlation function is a 
much easier object to calculate analytically as compared to the distribution. Such a 
calculation for ideal systems at zero absorption has actually already been performed 
in many interesting cases. In the particular case of a chaotic cavity an exact result for 
the correlation function (|33J) has been previously derived by us in ^nj|^]. Moreover, 
similar formulae can be derived for any model describing a one-particle quantum 
motion in a d-dimensional sample with a static disordered potential, see outline 
of the derivation in Appendix B. The only important physical assumption is that 
the disorder is such that all the relevant statistical properties of the system can be 
adequately described by the standard diffusive supersymmetric non-linear cr-model (or 
its lattice version). For a detailed discussion of the validity of this approximation (and 
its limitations) the interested reader is referred to the review j2H] and the book |74j . 

In the "zero-dimensional" case (chaotic cavity) the analytic continuation of (|33|l 
to complex O = ir can be represented generally as follows (see Appendix B for details): 

~ii 1 1 / d 2 d 2 \ ... . z' 2 + z" 2 + l . . 

^ «") = z n + {z n + 1)a + I [d^ + d^ ) n*) , *= 2z „ (37) 

where it is important that the function J-(x) depends on z' and z" only via the scaling 
variable x > 1. Its explicit form depends on the symmetry present (e.g. preserved or 
broken TRS), the following common structure being however generic: 



(38) 



T(K f 1 ^ f°HX f°H) /({A})e-^-^)/ 2 (x + A ) 

The real function /({A}) is the only symmetry dependent term here. In the crossover 
regime of gradually broken TRS it can be represented explicitly as follows [EI] : 

/({A}) = {(1 - Ag)(l + e- 2Y ) (A? - Ag)(l - e- 2Y ) 

+ 4y 2 K[(l Xl)e- 2Y + A|(l - e"^)]} (39) 

with 1Z = Ag + X'l + X 2 — 2AoAi A 2 — 1 and Y = y 2 (l — X 2 ), where y denotes a crossover 
driving parameter. Physically, y 2 ~ SE y /A is determined by the energy shift 5E y 



Scattering, reflection and impedance of chaotic waves 



14 



of energy levels due to a TRS breaking perturbation (e.g., weak external magnetic 
field in the case of quantum dots) (851 186) . Such an effect is conventionally modelled 
within the framework of RMT by means of the "Pandey-Mehta" Hamiltonian |87|. 
H = Hs+i(y /^/N)Ha, with Hs (Ha) being a random real symmetric (antisymmetric) 
matrix with independent Gaussian distributed entries. The limit y — > or y — » oo 
corresponds to fully preserved or broken TRS, respectively. 

Now we apply relation 1)36)1 to Eq. 1)37)1 and then perform the inverse FT to 
get V(u,v). Relegating all technical details to Appendix C, we emphasize here 
the most important points. The nontrivial contribution to the distribution comes 
from the second ( "connected" ) part of the correlation function 137)) whereas the first 
("disconnected") one is easily found to yield the singular contribution S(u)S(v — 1). A 
careful analysis shows that due to specific i-dependence given by Eq. (j3*S|) the above 
described procedure for the analytic continuation of the connected part of C(k, z") is 
equivalent to continuing J-(x) analytically and taking the jump at 

x = -x ± iO, with x = ±(u 2 + v 2 + 1) > 1 . (40) 

Thus, the nonzero imaginary part F(x) = IxilF(— x + iO) of the analytic continuation 
of ()38)l is determined at given x by the following integral 



l 



/({A})e-T( AlA2 - A »)/ 2 (a;_Ao) 



F{x) = JJ x °JL iMX * [(a?- m - 1) - (a,a: - (41) 

over the integration region B x = {(\\, A2) | 1 < \i,2 < 00, (AiA 2 — x) 2 < (\ 2 — 1)(A| — 
1)}, where the square root in ()38)l attains pure imaginary values. Taking now into 

account the following identity + F(x) = v~ 2 -^(x 2 — l)-^F(x), which is 

valid for x 2 7^ 1, we arrive finally at 

P(u,v) = j^- 2 Ux 2 - ^ l^Po(x) . (42) 

47r z ir ax ax 2ttv z 

This proves in general, cf. ()28)l . the statistical independence and uniform distribution 
of the scattering phase at perfect coupling. At arbitrary values of the crossover 
parameter y the obtained expression can be further treated only numerically. However, 
further analytical progress is possible in the limiting cases of pure symmetries and is 
considered below. We discuss explicitly the two most important cases of unitary and 
orthogonal symmetries. Clearly, the general scheme can be extended to the symplectic 
case with no principal difficulties. 



4- 1.3. Unitary symmetry. This case corresponds to considering the limit y — > 00. 
The nonzero contribution comes then from the second line in 1)39)1 where one can use 
effectively 4y 2 e~ 2y ( A 2 _1 ) — > <5(A 2 — 1) in the limit under discussion. This simplifies 
expression further to F(x) = /^dAi j\ (Al d _^ )5 fffi e" 7 ^ 1 ^/ 2 . Performing 
now analytical continuation ()40)) and making use of Im Ai _^. +i0 = —tt6(Xx — x). one 
gets readily F(x) = n j\d\ (x - A o ) _1 e - '''0 B-A( O/ 2 that yields the distribution [HI 

leading to equation ()31)l with (3 = 2 obtained earlier ^H] by a different method. 
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4- 1-4- Orthogonal symmetry. This case amounts to investigating (|39fl and (|41|l at 
y = 0. Fortunately, further simplifications are possible if one considers the integrated 
probability distribution W{x) = ~^{x 2 — 1)4-F(x) = dx P (x), which is a 
positive monotonically decaying function by definition. To this end, we note that it 
is useful to switch to the parametrization of |2] to carry out the threefold integration. 
The latter turns out to yield then a sum of decoupled terms and, after some algebra 
(see Appendix C), the result can be cast [22j in the following final form: 

x + 1 r i 

W{x) = — — fi{w)g 2 (w) + f 2 (w)gi(w) + hi(w)j 2 (w) + h 2 {w)ji(w) x (44) 

47T L J w— — ^ — 

with auxiliary functions defined as follows: 



Jt\t - w\ e-T*/ 2 r 

M " ,,= / ^ (l + t)3/2 

1 e~ 7 */ 2 

gt(w) = I dt- 



J\t - w\ p-^/ 2 
hi(w)= I dt Vlt [ 7+ (i- e -7) (7 t_ 2 )] 

Jl(U7) = / dt- 



\Jt\t-w\ VT+t' 

their counterpart with the index 2 being given by the same expression save for the 
integration region t G [0,w] instead of [w,oo). The interpolation expression 131( 1 
with = 1 compared to the above exact result shows systematic deviations in the 
nonperturbative regime of moderate absorption 7 ~ 1 (see |22j). Both formulae 
essentially coincide in the limiting cases of weak or strong absorption, when one can 
already use the more simple and physically transparent exact limiting statistics (|30|l . 

4-1.5. Impedance, reactance and the LDoS. Now we discuss statistics of the real 
and imaginary parts of the local Green function one by one. Let us consider first 
the distribution of the imaginary part v (LDoS). Having Pq(x) at our disposal, it is 
immediate to find the distribution function V v (v) of v by integrating out u in 128|) : 

V2 



V v (v) 



TTV 



3/2 







/ dqP 




Jo 





(45) 



This distribution is normalized to 1 and has the first moment unity at arbitrary 
absorption (indeed (v) = J Q dvvV(v) = 1 is automatically satisfied due to invariance 
of the integrand of with respect to the change «—►-). Explicit results for V v (v) 
as well as for V u {u) (see l(4"5|) below) can be readily obtained as Pq{x) at arbitrary 7 
has been already derived above. We refer the reader to the original publication |21j. 
focusing now on the physically interesting limiting cases of weak and strong absorption, 
where the behavior of the distributions is qualitatively different. 

At 7 <C 1, the exact limiting statistics l(3()ll for Po(x) can be used to perform the 
integration in Q45[). One arrives at the following result: 

a (l+f3)/2 v -(3+0)/2 e -a/4v ; y ^- a 

V v {v) cx { a 1 /2 w -3/2 j a<u<Cl/a (46) 

a (l+/J)/2 t ,-(3-/J)/2 e -a«/4 ) l/ a<U) 
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Figure 2. Distribution of the real u and imaginary v parts of the local Green 
function in the time-reversal chaotic cavity as absorption 7 grows. Dashed lines 
in (a) and (c) correspond to the limiting statistics 1491 and 1471 . respectively. 



where constant factors of the order unity are omitted. This result can be physically 
interpreted in the single- level approximation |88l I89| , when the main contribution to 
the LDoS v « v n — \^ n \ 2 ^{£n + jl 2 )^ 1 comes from a level closest to the given 
value of the energy (e„ — 2tt(E — E n )/A). Fluctuations of the wave function ip n arc 
mainly responsible for the exponential suppression of the distribution on the tails, 
while spectral fluctuations produce the intermediate bulk behavior with a "super- 
universal" |-power law, which actually does not depend on the symmetry at all. As 
7 increases, a number of levels contributing to v ~ ^ v n also grows ~ 7 and their 
contributions get less correlated. The resulting limiting form is almost a Gaussian: 



V v (v) 



exp 



— r W- 



1 v-' 

7^ 



a > 1 



(47) 



It shows a peak at v ~ 1 of the width oc 1/ y^y <C 1 , in agreement with (55] ■ 

Along the same lines, we may consider the distribution V u (u) of the real part u 
( "reactance" ) . One finds from l|28l) 



V u (u) 



1 



dqP Q 



<l 



(48) 



The limiting forms of V u (u) at weak and strong absorption follow readily as 

ti-^I + u 2 )- 1 , a<l 



V«(u)~{ . (49) 

^a/47re- Q ""/ 4 , a>l 

The Lorenzian distribution in the first line of Ij49(l is the consequence of the uniform 
distribution of the scattering phase [HO], as u ~ c °t f hi the limit considered. As 
absorption grows, one can see the crossover to a Gaussian distribution centered at 
zero, which is again a result of the central limit theorem. This type of behavior as 
well as a trend of V v (v) to the Gaussian distribution l|17)l. see figure 2, was recently 
observed in experiments on the cavity impedance |251 126| . 



4-1-6. Delay time at vanishing absorption and eigenfunction intensity. The Wigner's 
time-delay is one of the most important and frequently used characteristics of quantum 
scattering. It was recently realized (23 that the fundamental relation H28JI allows one 
to relate the distribution of the time-delay q in single-channel scattering (i.e. the 
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Wigner's time-delay) from a general lossless disordered system to the distribution 
of wave function intensities in the closed counterpart of such a system. Namely, 
consider the local eigenfunction intensity y = V\i/j n (r)\ 2 at a spatial point r of the 
closed system, with V denoting the volume of the sample and the index n numbering 
different eigenfunctions. Let V y (y) stand for the distribution of this intensity, where 
the statistics is sampled both over various realizations of the disorder in the system 
and over a certain small energy range around the point E in the spectrum, with A 
being the mean level spacing in that energy range. Define the dimensionless time delay 
t w = qA/2ir corresponding to wave scattering at the chosen energy E from a single 
channel attached to the point r of the sample. Denoting V w (t w ) the corresponding 
distribution in the regime of perfect coupling to the sample, Ossipov and Fyodorov 
managed to derive the following functional relation between the two distributions: 



Referring the interested reader for the details of derivation to we would like only 
to mention that at the starting point of the derivation the time delay is expressed via 
the reflection coefficient r in presence of absorption: t w = lim 7 ^o ^—*-> see @). Then 
one exploits in a clever way the scaling form of the distribution l|28|) remembering 
both the interpretation of the variable x in terms of r, see ()29f) . and the interpretation 
of v as the LDoS, the latter step providing a connection to eigenfunction statistics. 

On one hand, in the idealized situation of zero absorption statistics of delay times 
of all sorts were studied intensively in the framework of the RMT approach, various 
exact analytical results being available, see |ini E??ll^l lMll?mi?^ll?^llMll??KllM | 
and references therein. Those were successfully verified in numerical simulations of 
chaotic systems of quite a diverse nature, see jH3 USD EH)- Since phase shifts and 
delay times are experimentally measurable quantities, especially in a single-channel 
reflection experiment |U HOI ITUU1 ITU1 ITU21 ITU5| . the relation opens a new 
possibility for experimental study of eigenfunctions. 

On the other hand, one can use the existent knowledge on eigenfunction statistics 
|1U41 129) to provide via (|50|) explicit expressions for time-delay distributions. In this 
way one can e.g. recover those for chaotic systems of all symmetry classes obtained 
previously by diverse methods in various regimes of interest gUl EU Ej . Of particular 
interest is the predicted multifractal behaviour of the negative moments of time-delays 
in the vicinity of the Anderson localisation transition [M5| : 



where L stands for the system size at criticality, and D q are anomalous (multifractal) 
dimensions of the eigenfunctions. Such a behaviour was indeed discovered recently 
|1U5) in numerical simulations of the disordered lattice Hamiltonians at criticality. 

4- 1.7. Anderson transition as phenomenon of spontaneous breakdown of S-matrix 
unitarity. Actually, absorption in disordered systems plays not only purely technical, 
but also a conceptually important role in revealing the mechanisms behind the 
Anderson localisation transition. Lets us shortly discuss a possible qualitative 
behaviour of the PNR r = 1 — r in a scattering system formed by a single perfect 
channel attached to a d-dimensional disordered sample at the vicinity of the point of 
the Anderson delocalisation transition /i c (the mobility edge). Here we denote by /x 
an effective parameter which controls the transition in the infinite sample, with states 
being localised (extended) for fi > [i c (respectively, jU < pL c ). 



(50) 




cx L~ qDq+1 



(51) 



Scattering, reflection and impedance of chaotic waves 



18 



For a sample of finite size L the PNR is a function of three parameters: absorption 
r, size L, and disorder strength fi. In the insulating phase fi > y, c the system is 
characterized by a localisation length £ which diverges when fj, approaches the critical 
value. A natural scale for the absorption is played by a parameter 8 oc T£ d , i.e. the 
ratio of the imaginary shift in energy to the mean level spacing for sample of a typical 
size L oc £ (localisation volume) . 

Consider the weak absorption limit 8 — > 0. It is clear that the wave incoming 
through the incident channel effectively explores only a part of the sample of the order 
of localisation volume £ d , being exponentially small elsewhere. Under this condition 
it is clear that the two limits: T — ► and L — > oo should actually commute and can 
be taken in arbitrary order. We know that in the limit r — > the PNR behaves as 
r ~ Tq — ► 0, see (JSJ. Moreover, exploiting the relation (| 5 C)|) . and remembering that 
(\ip n (r)\ 2k ) oc it is easy to see that all negative integer PNR moments should 

behave in the infinite volume limit L — > oo as (r~ k ) ~ p — *£ — . With a little bit more 
work one can suggest a qualitative picture for the PNR probability distribution in the 
localised phase. Namely, the distribution should have a powerlaw tail V(t) oc 8/t 2 
extending through a parametrically large domain 5 < r < 1, and should decay very 
fast towards zero for both r <C 5 and r — > 1. When absorption vanishes 5 oc T — > 
such a distribution collapses to the Dirac 8— function: V(t) oc <5(r), but in a very 
nonuniform way. We may conclude that in the limit of vanishing absorption 5-matrix 
unitarity is indeed recovered, and in this sense we can associate the localised phase 
with the phase of unbroken symmetry. 

In contrast, in the delocalized phase the incoming wave explores the whole sample 
volume. It is natural to think that whatever small (but fixed) is an absorption rate 
r, in the limit L — > oo a finite portion of the incoming flux will be absorbed in the 
sample and will never come back to the incident channel. In particular, we may 
expect limr^o bniL^oo r(r, L, /i) = Too(^) > as long as fi < fi c . From this point 
of view the Anderson transition acquires a natural interpretation as the phenomenon 
of spontaneous breakdown of 5-matrix unitarity. Informal arguments in favour of 
such a behaviour are based on a picture of the transition as described in terms of a 
functional order parameter developed in detail in pnj[S21, see also earlier results in [71) 
and [106 . Namely, the distribution function V(r) is expected to remain a non-trivial 
finite- width distribution even when T — > 0, provided the latter limit is taken after the 
infinite volume limit L — + oo. Clearly, more work is needed to substantiate this claim, 
as well as to clarify critical behaviour of Too(/i) as long as /i — > /i c . 

Note finally, that if the limit r — > is taken first, then for /i < fi c PNR in large 
but finite sample should again scale with the system size L as (r~ fc ) ~ C(fi)T~ k L~ dk , 
cf. H51|l. The coefficient C(/z) is however expected to diverge when fj, — » fi c . This 
divergence should be related to the properties of eigenfunctions via equation H50J1 • 



4-1.8. Arbitrary coupling to the channel. The general case of arbitrary transmission 
coefficient, T < 1, can be mapped |1()71 11081 1951 ITT)) to that of perfect coupling by 
making use of the following relation 



S T=1 = S-VT^T 

between the corresponding scattering matrices. Equation (|52|l is known from the 
Poisson kernel theory |1()7| . Due to an additional interference between incoming 
and directly back-scattered waves, the scattering phase 9 acquires a non-uniform 
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distribution and acquires statistical correlations with x (or r). However, the joint 
probability density P(x, 9) is again determined by the function P as follows [21] : 

P(x, 6) = ^P (xg - V(* 2 ~ 1)G? 2 - 1) cos 9) , g=^-l>l (53) 

This relation can be obtained by a straightforward evaluation in the parametrization 
<|52[1 of the corresponding Jacobian. The integration over x immediately yields 
the scattering phase distribution. In the particular case of vanishing absorption 
r — ► 1 and Pq(x) — > x _2 (5(l/a;). This readily gives P(x,9) — > p(9)Pq(x), with 
the phase density 2np(9) = (g — \J g 2 — lcosfl)^ 1 found earlier [HHl (see j^H for 
the corresponding numerical study). As another example we consider the reflection 
coefficient distribution in terms of Pq(x) is given at arbitrary T by |17|: 

Distributions of r and in microwave cavities were recently studied for different 
realizations of 7 and T in |l()j . the excellent agreement with the theory being found. 



z{g — y g — i^rcosi 



1 - r 



(54) 



4-2. Beyond single channel 

4-2.1. Reflection coefficients and PNRs. The starting point of our analysis in this 
section is the following convenient representation |17j for the diagonal elements S CC (E) 
of the scattering matrix, cf. (JTJ : 

S CC {E) = \~w°!{E) ' Kc[E) = ■ ^ 

where H^. s — H — ^ V b V b1i is now the c— dependent non-Hermitian operator. 

In view of H cS = H - \VV^ = H c cS - {V C V C ^ one can treat V C V C ^ as a rank one 
perturbation with respect to H-eff- ^ n this case the following general relationship 
(Dyson's equation) is valid for the corresponding resolvents |38| : 

' ' ' ' V l—y* 1 (56) 



E-H cS E-H c cff 2E-H c cS l + \K c {E) E-H c eS 

which can be proved by expanding the l.h.s. in a power series with respect to 
(E — TL^ S )~ 1 V C V C ^ and summing then up the resulting geometric series in K c . 
Substituting this equation in S cc = 1 - iV ct (E - H cff )~ 1 V" c one gets S cc = 
1 — 2iK c — 2K c (l + \K C )~ 1 K C which is equivalent to (|55|l . It is also worth noting 
that a representation like (|55|) is valid for an arbitrary mxm submatrix standing 
on the main diagonal of S with obvious replacement V c by the Nxm matrix 
y(m) _ (y ci , . . . , V Cm ), and corresponding changes for K c and Ti^ s . In the particular 
case m = M, i.e. the full 5-matrix, one recovers from (0) and i|56|) . 

The relation l|55Jl reduces the problem of evaluating the statistics of S cc , and hence 
the reflection coefficient / PNR in a given channel c to calculating the joint probability 
density V{u c , v c ) of u c = KeK c and v c = —lmK c , with K c standing for the particular 
diagonal entry of the resolvent of the effective non-Hermitian Hamiltonian 7ijj ff from 
15511 . Moreover, a uniform absorption within the sample can again be taken into 
account by a purely imaginary shift of the scattering energy E — ► E + = E y> 
all further steps being fully in parallel to those of section 4.1.2. The most pleasant 
feature of this approach is that it is very straightforward to include open channels in the 
derivation of V(u c ,v c ) within the supersymmetry method. All important properties 
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of this distribution, in particular, the relations Q2H|). <|41|) and (|42|l retain its validity, 
with x being naturally replaced by x c = ^r(Uc + v c + !)■ The corresponding function 
F(x c ) then follows from (|41|l by multiplying there the integrand with the "channel 
factor" (which originates from the imaginary part of Ti^ff ) 

9b = Y b ~ 1 ( 5? ) 

which accounts for arbitrary coupling to all other channels save for the given perfectly 
open one c. Arbitrary coupling for the remaining channel c can be considered by means 
of i|53|) . providing us finally with the general distribution of the reflection amplitude 
and phase. In this way one can also obtain explicit distributions for many interesting 
situations [T7 , including the cases when the effects of Anderson localisation start to 
play an important role. In particular, in the case of broken TRS one can calculate 
the PNR distribution for a single channel attached to an edge of a piece of quasi- 
one-dimcnsional disordered medium of a given length L, with the opposite edge being 
either closed, or in contact with a perfectly conducting multichannel lead. 



n,({A» = n 



(_9b + Ap) 2 

+ A 1 A 2 )-(A?-l)(Af-l). 



4-2.2. Reflection eigenvalues and thermal emission. The exact result for the 
distribution function P(r) of reflection eigenvalues valid for any number of arbitrary 
open channels and arbitrary absorption is known only for the (3 = 2 case, being recently 
calculated by two of us [TS], generalizing earlier perturbative results (known for all 
(3). That uses the method developed for studying the proper time-delay distribution 
96 in the ideal lossless system. The later distribution considered at finite absorption 
by means of © has a sharp cutoff at q — T _1 due to a finite value of P(r) at r = 0. 
This fact makes the interpretation of q as delay times at strong absorption questionable 
|16j . since intuitively one expects a generic exponential suppression at large values of 
delay times q ^> Indeed, for the time St a wave-packet oscillating in the cavity 

with a frequency A/2n on average experiences (A/2%)6t collisions with the walls, 
yielding the probability T < f,(A/2w)6t to be absorbed into one of parasitic channels 
(T<£ <C 1 being the transmission coefficient). The total reflection is then estimated as 
R ~ (1 — T < f } (A/2n)5t) M <£ , giving e~~ rSt in the absorption limit of the fixed absorption 
width r = M ( f,T^ ) A/2'K as — > oo and — * 0. It is natural, therefore, to define 
alternatively the positive definite matrix Qr of reflection time-delays as follows |16| : 

q r = -T- 1 in r 7 = -r- 1 in(i - rg 7 ) . (58) 

One finds easily the connection Vn(q r ) = e~ Vqr, P[j;{]- — e~ r<?r )] between the 
corresponding distributions T , R(q r ) and V(q) of proper delay times (eigenvalues of Qr 
and Q 7 , respectively). Both Qr and Q 1 reduce to the same Wigner-Smith matrix Q 
in the limit of vanishing absorption. The difference between them becomes noticeable 
at finite T. Still both distributions coincide up to the time <C They start to 

differ at larger times, when V(q) has the cutoff whereas pR{q 3> T _1 ) cx e~ Tq is 
exponentially suppressed. 

The exact result for the f3 = 1 case at arbitrary M is still outstanding. In the limit 
of the large number of equivalent channels, M — ► oo, an exact result can be found at 
arbitrary absorption and coupling |18| . and the perfect coupling case T = 1 has been 
known for some time |109M110"| . In contrast to the few-channel case, when any value 
< r < 1 is permitted, the distribution density in the present case is non-vanishing 
only in a range < r m ; n < r < r max < 1, and is the same for all f3. Referring the 
reader to JT& for explicit results, we mention their application for thermal emission 
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from random media. Registration of n photons in the frequency window 5Q during 
the large time t 3> 1/5VL yields the negative-bimodal distribution of photocounts with 
v = Mt5n/2w degrees of freedom; see |111| . In his seminal paper Beenakker |109) has 
shown that the quantum optical problem of the photon statistics can be reduced to 
a computation of the S matrix of the classical wave equation. In particular, chaotic 
radiation may be characterized by the effective number v c g degrees of freedom as 
follows |109| : v e ff_/v = (1 — (r)) 2 / ((1 — r) 2 ) < 1, with v c s = v for blackbody radiation. 
This ratio is given by v e s/v = (p/ s + T) 2 /[j 2 + 2(j s + T)] at arbitrary transmission T 
and absorption 7 S = 2ttT/AM, implying thus that saturation to the blackbody limit 
slows down, ^ « 1-2 (1 - T) at 7 S > 1, for transmission T < 1 [JSj. Finally, 
due to a duality relation |1U9I I112j of a linear amplifying system to a dual absorbing 
one (related to it by inverting the sign of T) there exists a further link of the analysis 
presented here to the rapidly developing field of random lasers |109H110llll2llTT5llll4| . 



4-2.3. Off- diagonal entries of the Green function. As we have seen, the statistics 
of real and imaginary parts of diagonal elements of the Green function can be very 
efficiently studied in the framework of the supersymmetric approach and may have 
various physical interpretations. The off-diagonal entries Gij(E^) — (i\(Ey— Ti. c s)~ 1 \j) 
are of considerable importance as well. It can be easily understood that W — \Gij\ 2 
is essentially the wave power transmitted from a source at site i inside a random 
medium to a receiver at site j. Statistics of such an object is much more difficult to 
study in general. Presently the most studied case |19| is /3 = 2 symmetry class under 
an additional assumption that both receiver and source are very weakly coupled to the 
medium as to ensure those couplings do not contribute essentially to the resonance 
broadening. This means the broadening is induced purely due to losses elsewhere 
in the medium. Technically the latter requirement amounts to vanishing coupling 
amplitudes Vf and Vf for all damping channels c = 1 . . . M. Assuming further the 
RMT statistics for Re7i c ff, it is easy to see that the damping matrix Im7i c ff can be 
chosen diagonal and simply such that both entries (VV')u and (VV')jj are vanishing. 
For such a model the distribution of the scaled transmitted power w = (-§-) 2 W can be 
found explicitly as: 



v(w)= ( A + 

aw aw 



M 



sf\ + w + X 7A j-r g c + A 



./i j- „, _ \ ii 



2yl + w 7-1 Vi + w - A g c + Vi + w 



(59) 



In fact, the remaining integration can be performed for two interesting cases: (i) all 
equivalent dissipation channels g c — g in the absence of uniform losses 7 = and 
(ii) no internal channels of dissipation g c — > oo in the presence of uniform absorption 
7 > 0. Here we present the formula only for the latter case |19| : 

V(w) = 1 — 5 1 [ 1 l( w + 2)-(w-2-2( 1 l/ 7 )coth 7 )(l + lw )] (60) 



with the shorthand j w = y\fl + w. The distribution of the transmitted power for 
the case of absorptive media with preserved time-reversal invariance is not yet known. 
First two moments of that quantity were calculated recently in |20| . 

Let us finally mention that a closely related question of statistics of intensity of 
waves emitted from a permanently radiating source embedded in a random medium 
was the subject of many studies in recent years. Some useful references can be found 
in Section 7 of [2H], see also the relevant review |114| on random lasing. 
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5. Conclusions and open problems 

For wave scattering in open chaotic and/or disordered systems with uniform losses, 
we have discussed various statistics on the level of both correlation and distribution 
functions. The overall exponential decay due to uniform absorption is the generic 
feature of any correlation function reduced to a two-point spectral (resolvent) 
correlation function, that follows simply from analytic properties of the latter in 
the complex energy plane. For fully chaotic systems with or without TRS, we have 
calculated exactly energy correlation functions of complex impedance and 5-matrix 
elements at arbitrary absorption and coupling. The corresponding enhancement 
factors have been also discussed in detail. The result for 5-matrix correlations in 
the case of broken TRS completes the well known one of preserved TRS. 

To study distribution functions, we have described the novel approach to the 
problem by deriving a kind of a fluctuation dissipation relation H36|l of quite a general 
nature, which relates distributions at finite absorption to arbitrary order correlations 
at zero absorption via a nontrivial analytic continuation procedure. Correlations can 
be efficiently studied in the framework of the supersymmctric nonlinear cr-model. 
In the zero dimensional number of explicit results is provided for quantities 

characterizing open chaotic systems both from "inside" (LDoS, the complex impedance 
and the (local) Green's function) and from "outside" (PNRs, reflection coefficients and 
time delays). This cr-model mapping provides an attractive possibility to include into 
consideration open disordered absorptive systems beyond the usual RMT treatment. 

Finally, let us briefly discuss an (incomplete) list of problems deserving, from our 
personal point of view, further investigations both for scattering systems with and 
without absorption. Within the domain of RMT applicability, the most outstanding 
problems are of course those requiring evaluation of four-point correlation functions of 
the resolvents, e.g., partial scattering cross-sections. These are still not known even for 
the simplest case of broken TRS. Another challenge is to find a probability distribution 
for the multichannel 5-matrix at finite absorption, thus nontrivially generalizing the 
Poisson kernel j9QJ. More work is required to understand an interplay between the 
statistics of complex resonances and the corresponding bi-orthogonal eigenfunctions 
[381 169| . the question being of particular relevance for lasing from random and/or 
chaotic media |115II116| . Spatial characteristics of internal parts b c = (E — 7i e fi)~ X V° 
of the scattering wave function are very interesting by their own and are closely related 
to fluctuations in transmitted power discussed above in 4.2.3. 

Apart from that, last decades it was realized |117j that it may be required to 
consider other symmetry classes (beyond the three of Dyson) which are relevant, e.g., 
for systems involving superconducting elements. The corresponding scattering theory 
is reviewed in |118j . To reconsider many problems discussed in the present review 
taking those (and other related) symmetries into account should be both interesting 
and informative (see |24j as a recent example). 

At last but not least: understanding scattering in disordered systems beyond the 
universal RMT regime, taking into account, in particular, the effects of Anderson 
localisation, is a rather promising area calling for more systematic research. 
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Appendix A. Statistics of the local Green function in the GSE 

We start with expressing the dimensionless LDoS (the imaginary part of the local 
Green's function G{E + ^T; r; r) in units of A) in terms of the eigenfunctions and 
eigenvalues of the Hamiltonian H with underlying simplectic symmetry as: 

v{r) = -^I mG (i? + ir;r,r) = [|^„(r)| 2 + |0„(r)| 2 ] (jr _J,. lp9 . (A.1) 

Here tpn( r )> 0n( r ) stand for the local amplitudes of the two eigenfunctions 
corresponding to the (double-degenerate) eigenvalue E n . It is convenient to consider 
scaled eigenvalues e n — irE n /A, defining absorption rj = 71T/A = ^7. Replacing H 
with 2N x 2N random GSE matrix, we follow the idea first suggested in |119| and in the 
first step exploit the well-known fact that in the limit N 3> 1 eigenvector components 
i/j n (r), 0„(r) for different values of n and r can be treated as independent, identically 
distributed complex Gaussian variables. The corresponding joint probability density 
can be written symbolically as V n (ip; tp)V(ip; <f>) = e~ N ^ ^V(ip;(f>) where we 
introduced a shorthand notation tp^ip — 52 n=1 |^„(r)| 2 (and similarly for <jy<j>), with 
the measure T>(ip; 0) being understood as the (normalized) product of differentials of 
independent variables. 

The Laplace transform F(s) — / °° e~ ps V p (p)dp of the probability distribution 
function V p for the normalized variable p = Nirv(r) can be then written as: 

PM / [ t> ( 1 M\t>{( 1 Ax\\T\ ^(IV'»(r)| 2 + |^(r)| 2 ) \ 

F{s) = { / V n ty;$)V(M>;$)) I I cxp K - '- ) A.2 

where the angular brackets stand for the averaging over the joint probability density 
of all eigenvalues e n and the limit N — > 00 is to finally taken. After performing the 
Gaussian integrals, we therefore arrive to the following representation: 

n.)-**(fl K g - e "> 2+ " 2 ] 2 \ (A.3, 

"-~\,ii[(<-<..) s + ') 2 + 2»»]7 l «„ l 

Introducing the product 7Tat(A) = II^Li ~ e «) we see that the characteristic 
polynomial of the GSE matrix H is simply 7r^(A). The formula i|A.3|) can be 
conveniently rewritten in terms of this polynomial as: 

F{s)=^ / ^i±j^izM \ , M = V^T^ (A.4) 
\ < /2 (c + iM)< 7 (e - V) / {en} 

with (3 = 4 for GSE. In fact, in such a form the formula retains its validity for GUE 
(3 = 2 and GOE /3 = 1 cases as well. 

The problem of evaluating ensemble averages of products and ratios of 
characteristic polynomials of random matrices attracted a lot of research interest 
recently, both in physical and mathematical community. For /3 = 2, the general 
problem was solved in |12(J1 I121| , and the most complete set of formulae available 
presently for /3 = 1, 4 can be found in the recent paper by Borodin and Strahov |122j . 
When addressing the most interesting case j3 = 1 they were, unfortunately, able to 
consider only integer powers of characteristic polynomials, whereas (IA.4II obviously 
requires knowledge of half-integer powers. In this sense the corresponding random 
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matrix problem for (3 = 1 case is still outstanding. Borodin and Strahov derived the 
following explicit expression for the (3 — 4 average featuring in (|A.4|I : 



lim (%S^4) nu^-/^) Pfd (A . 5) 

JV^oo \ T7 2 N (p 1 )7r 2 N (f3 2 ) I GSE (ai - a 2 )((3i - (3 2 ) 

assuming that Im/3i > 0, Imp2 < 0. Here, the matrix C has the following structure 

C=(S I),withi=(_° a * ),D=(° d ^) andPf 

stands for the corresponding Pfaffian, Pf(7 = V det C. The entries of the matrix C 
are given explicitly by: 

a = - ^ T sin (ai - a 2 )t, d = 2.1— — j (A.6) 

e i(/3i-Qi) e i(ai-/3 2 ) gi(/3i-a2) e i(a 2 -/3 2 ) 

&n = — a , bi2 = — , b 2 i = — , b 22 = — (A.7) 

Pi — cti ai — p 2 pi — a 2 a 2 — p 2 



Restricting our consideration in l|A.4|) for simplicity to e = E = we have therefore 
ai = in, a 2 — —in,(3x — i\i,(3 2 — —ifx. Substituting this to (|A.6fl — (|A.7f> and then to 
l|A.5|) , we find after straightforward computations the Laplace transform of the LDoS 
probability density: 

Fgse(s) = F GUE (s) + SF(s), (A.8) 
FgUe{s) = TJi ( e ~ 2( "~' ,) ('7 + M) 2 " e-^Hn - M ) 2 ) (A.9) 

SF(s) = ^(v 2 ~ M 2 ) 2 e- 2 " (l + ^) jT 1 j sinh(2r,t), (A.10) 

We note that (i) the only s— dependence in the above expression comes from fi = 
\frj 2 + 2ns and (ii) the formula (|A.9J| concides with the Laplace transform of the 
LDoS distributuion for (3 = 2 case, see e.g. |88l 1119) . 

The remaining job is to find the part SF(p) of the LDoS probability density 
corresponding to inversion of the Laplace transform in l|A.10f) . which gives 

SF(p) = ^1 Li/2 e -n( P /2+2/p)l x 1 [X f 1 ± sinh (2nt) (A.ll) 
dp L J 4 V irn Jq t 

Knowledge of this function allows immediate restoration of the corresponding 
probability density Il-Sl'l) for the main scaling variable x due to relations discussed 
in the main body of the present paper (recall that n — "f/2). 



Appendix B. Disordered systems of arbitrary dimension: 
nonlinear er-model derivation of expressions (|37|) 

Let H be a (self- adjoint) Hamiltonian describing a one-particle quantum motion in a 
d-dimcnsional static disordered potential. Among microscopic model Hamiltonians 
ensuring validity of the non-linear sigma-model description of such a system the 
simplest choice seems to be the Wegner's N— orbital model [1231 1124] . or its variant 
due to Pruisken and Schafer |125j . Physically the models are equivalent to the so- 
called "granulated metal" model [23]. One can visualize it by considering a lattice of 
L d sites (d standing for the spatial dimension of the sample) , each site being occupied 
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with a metallic "granula" . The motion of a quantum particle inside each granula is 
assumed to be fully " ergodic" , and as such the Hamiltonian of the individual granula 
should be adequately modelled by a random NxN matrix H of appropriate symmetry, 
provided we consider the limit N — > oo. The quantum particle is also allowed to tunnel 
between the neighbouring granulae, the process ensuring a possibility of nontrivial 
diffusive motion along the lattice. Thus, the Hamiltonian H of the system as a whole 
has a form of large matrix of the size N x L d , consisting of coupled matrix blocks 
Hk, k — l,2,...,L d of the size NxN. For example, for the simplest quasi-one 
dimensional sample d — 1 such a matrix will be "block-three-diagonal" . 

Being interested in scattering, we should provide a way to incorporate description 
of external leads (or waveguides) attached to the sample. In the framework of the 
present model attaching an M^-channel lead to the block at site k is done by replacing 
the corresponding "intragranula" matrix block Hk, with its non-selfadjoint counterpart 
H k - if fe , where NxN matrix f fc = 2 diag(/4*° , . . . , k { ^\ , 0, . . . , 0) > 0, and M k 
staying finite when N — > oo. Assuming the absorption width T to be uniform and 
identical in all the granulae, we are interested in deriving the joint probability density 
V(u,v) of real u and imaginary v parts of the local Green's function in a state 
belonging to granule at the lattice site I: 

G(jWJ l \E 7 ) = (i«|(£? 7 -H+lf)- 1 ^)) . (B.l) 

Following the general strategy (see Section 4.1.2) we recover V(u,v) from 
the correlation function C(z',z"), see (|34H . The calculation of the corresponding 
correlation function is a straightforward extension of the "zero-dimensional" procedure 
employed in 0O| for (3 = 2 and then in jH^ for f3 = 1 (as well as for the whole 
crossover) to the present d-dimensional situation and will not be repeated here. The 
emerging supersymmetric nonlinear er-model on the lattice is described in terms of 
the supermatrices Qk, k = 1, . . . L d parametrized as Qk = — iTfeAT^T 1 and interacting 
according to the "action" (see e.g. |%T]) 

S{Q} = -y E StrQ^ + I^^StrOkA (B.2) 

<k.j> k 

where the first sum goes over pairs < k, j > of nearest neighbours on the underlying 
lattice, A stands for the mean level spacing and t stands for the effective inter-granule 
coupling constant, which is the main control parameter of the emerging theory. We 
start with outlining the procedure for (3 = 2 symmetry when Q-matrices involved in 
the calculation have the smallest size 4x4. Assuming for simplicity that the granule 
/ does not have a channel attached to it directly, the correlation function C(z',z"), 
see and of resolvents of the (scaled with Green's functions (|B.1() is 

expressed in terms of the Q-integrals as follows 

L d 

Cp=2{z',z") = IJ / MQk)^l{Q k )e- s ^^ T SAet- l {Uj 1 +Q l )\ Ji=j2=0 (B.3) 
k=i J 

where the supermatrix U depends on the variables z' , z" and sources J\, J2 via 

U = \{z' - iz")(l + A) + \{z' + iz")(l - A) + diag(0, J 1} 0, J 2 ) . (B.4) 

The "channel factor" J 7 ^ (Qk) = Tic^i Sdet -1 (1 + i«c fc " ) Q fe A) originates from coupling 
to continuum, where k^) — Q if no external channels is attached to a given granule. 
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The factors Tm appearing explicitly in Section 3.2 as well as in l|57|l are just the zero- 
dimensional analogue of Tj^ . Following |81l l%2~] , we find it convenient to introduce 
the function 

Ld r 

Y m (Qi)=H /d/x(Q fc )^(Q fe )e- s «> . (B.5) 

k^l J 

Employing the standard Efetov-type parametrization of the matrices Q, one finds that 
actually Y">(Q) must be a function of only two commuting variables: 1 < Ai < oo 
and — 1 < Ao < 1. Then integrating out all the remaining degrees of freedom in quite 
a standard way (see 0U| f° r more detail), we arrive at the representation l)37|) where 

^=2(z',z") = fV f dA ^ + A y(0(A Aq) (B6) 
h J-i (Ai - XaV x + Xi 

the main scaling variable of our theory x = -^priz 12 + z" 2 + 1) being introduced. 
This immediately verifies the structure of the important formula (|37|) for a general 
d-dimensional (3 = 2 system with attached scattering channels. Of course, ability to 
work out explicit expressions for the probability density V(u, v) crucially depends on 
the availability of the function Fw(Ai,A). This function has a very simple form in 
" zero dimension" , physically equivalent to a system consisting of a single granula: 

M . 

rW(A 1 ,A ) = e -^ Al - A °)/ 2 n^ ± T £ (B.7) 

Li 9- + A i 

where we have assumed that the total number of external channels attached to the 
system is M. Each channel is characterised by its own effective coupling constant 
g c = ^(k c + k^ 1 ) > 1, the transmission coefficient being T c = 2 +1 . Weak-localisation 
corrections to zero dimensional results can be in principle found systematically 
following the procedures of |1(J4) . and non-perturbative results are mainly available in 
quasi-lD systems, see e.g. [126 , in the limit of weak absorption, 7 <C 1. 

Treating systems of /3 = 1 symmetry class follows exactly the same steps as 
outlined before, although in this case the supermatrices Q are of the size 8x8 with 
the corresponding doubling the dimension of (|B.4J| . Following 84, we consider the 
whole crossover between (3=1 and (3 = 2 symmetry classes. Due to the TRS breaking 
perturbation controlled by the parameter y (Hk is no longer symmetric but still 
Hermitian), the action l|B.2(l (with ,9=1) acquires additionally the symmetry breaking 
term S y {Q} = \y 2 J^k ^ T ( T 3Qk) 2 , with T3 being the Pauli matrix. As usual, explicit 
integrations are much more cumbersome and require more work, some very useful and 
helpful relations can be found in |86j. Exploiting the parametrization suggested there, 
the function Y^>{Q) turns out again to be dependent on three variables 1 < Ai,2 < 00 
and — 1 < Ao < 1, and the analogue of l|B.6(l reads 

T p=1 {z\z")= d\J dX 2 ^ + = y(')(A 1 ,A 2 ,A ) (B.8) 



where K = A§ + X\ + X\ - 2A AiA 2 - 1 and \x Y a = AiA 2 ± \J {X\ - 1)(A| - 1), and in 
zero dimension: 

M . 

Y^iXu A 2 , Ao) = f({X})e-^^>y2 -Q + *° _ . (B.9) 

c= i V(9c +Mi)(5c + M2) 
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Here, /({A}) coming solely due to the e~ Sy ^^ term is explicitly given by i|39|) . 
Integration variables /2i 2 make a connection between the "radial" parametrisation 
of Q of [72] and the "angular" parametrisation of [2] ■ This will be utilised in the next 
Appendix to get explicit result for the zero-dimensional case for (3 = 1. 

Appendix C. Analytic continuation (I40H and GOE result (I44[) 

We prove here equivalence of the analytic continuation of the "connected" part of 
C(z', z") in z" -> -t> ± iO JSBJl to that of ^(x) in x -> -x ± iO (gOJ). Let us consider 
in detail first the simplest case of (3 = 2 symmetry class in zero dimension. Then 

1 1 / d 2 d 2 



dA i / 1 — la _i_ ZJJ2 , <wm — TT C^- 1 ) 



1 



C(Z ' Z ) = z' 2 + (z" + l) 2 4\d 

rl dA 2z"e-y( Al - A °)/ 2 
_ x Ai - A z' 2 + z" 2 + 2z"Ai + 1 

where we have kept in the integrand only the term remaining nonzero under the 
differentiation. The FT C(k, z") with respect to the first argument reads as follows: 

7r e -l fc IO" +1 ) 1 / d 2 \ 

7_i Ax - A V^" 2 + 2z"Ai + 1 
Now we continue this result analytically in z" and calculate the jump (|36|l on the 
negative real axis z" — > — w ± iO, with u > 0. One gets readily 

2 dAl /' dA cos(fc^^^^)^ e _ 7(Ai _ Ao)/2 
1+id y_! Ai - A V2wAi — 1 — v 2 
Performing now the inverse FT one arrives at 

a dAi t -^-e-^- x ^/ 2 S(u 2 + v 2 + l- 2vXi) 2v . (C.4) 

Making now use of S(u 2 +v 2 + l-2vXi) 2v = ■^Q(2\ 1 v-u 2 -v 2 -1) = <S(Ai™x), with 
x = ^{u 2 + v 2 + 1) from (JSSJl, one can immediately integrate (|C.4|I over Ai, yielding 

fx—i t^ 1 e^ lt / 2 for the second line, that is in exact agreement with expression (|43|l 
following from the analytic continuation l|4(J|) . 

The proof for the (3=1 case proceeds along the same lines, explicit expressions 
being however more lengthy. We omit it, considering instead the less trivial derivation 
of expression l|44|) . It is important to stress that for (3=1 the nonzero F(x) (|41|) at 
given x > 1 is determined by the integration region B x = {(j5i, /S2) | x < jli < 00, < 
A 2 < x}, rather than by a single point Ai = x as in the (3 = 2 case. It is convenient, 
therefore, to choose £ix,2 from (|fj.8l) as new integration variables. Actually, they are 
related to those from |2] (which are already used in (|16|) ) as follows: 



Ml,2 = 



3(Ai,a-l) = | 



AiA 2 ± ^(A 2 - 1)(A 2 - 1) 



/i = |(l-Ao). (C.5) 
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with the pre-exponential factor in (| 1 fc>|) being the corresponding integration measure. 
Expression l|41|) at y = acquires then the following form (with w = ^i^): 



o Jw Jo + Mi)M2(l + A*2)(mo + Mi) 2 (Mo + M2) 2 

xc -7(mi+M2+2mo)/2 W + ^q 

- - M2) 







- 1)^2(^2 - l)(^o + ^i) 2 (^o + ^2) 



2 



^/(l + wi)(l + w 2 ) 

where the second equality comes after the scaling /xo,i,2 = w ^0,1, 2- It is convenient 
to consider now the function 4-F(x). An essential observation is that this function 
contains no contribution coming from the derivative of (|C.6(1 on the upper limit due 
to vanishing of the integrand at v a — 0, — . We represent it in the following form: 

= ( 1/w Avo r dui r ^ 2 e-^^x {W2 ) (c 7) 

dx J Ji J Q yjvx{vx ~ \)v 2 (vi - 1)(1 + wi)(l + wv 2 ) 

where a rational function X contains all other terms resulting from the differentiation. 
It is a crucial fact that X can be further decomposed into partial fractions with respect 
to vq, yielding 

0^(7 + ^12) . ih + Cidu 



x{v Q ,v l ,v 2 )= ^(-ir +i 



(vq + Vi) 2 Uq 



(C.8) 



with functions a, = Vi(i>i — 1)(1 + wvi), b% = 1 — 2(1 — w)vi — 3wvf, Cj = 1 — 2i>i — wvf 
for i = 1, 2 and di2 = [(1 + wi)(l + wz/ 2 )] _1 . Substituting (|C.8|I in (|C.7|I . one readily 
sees that integrals over i>o can be easily performed while remaining integrations over 
^1,2 get completely decoupled in each term of the sum, leading finally to 1441) . 
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